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inflaton mass much smaller than the Planck mass. 
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I. INTRODUCTION 


Cosmology is a well understood branch of physics today. The data gathered since the early 
1990s has allowed us to study thoroughly the origin of the observable Universe. One of the 
most important theoretical achievements is the inflationary scenario^. Cosmological inflation 
was introduced by Starobinski [3] and Guth [1] in the early 1980s in order to explain why 
the observable Universe is so homogeneous and flat. It has been known since 1929 [5] that 
the Universe is expanding, but if we trace this expansion backward in time, we see that 
the present-day sky was made of many causally disconnected patches. This either means 
that the homogeneity we observe today cannot have been produced by any interaction in 
the past, or, one can argue that the Universe was actually smaller than we expect from our 
simple extrapolation and it therefore nnderwent some kind of accelerated expansion in a very 
short time just after its creation. This brief expansion is called cosmological inflation, and is 
generally assumed to have occurred from about 10“^® to 10“^^ seconds after the (presnmed) 
origin of the Universe. 

The large amount of data collected by experiments such as COBE, BOOMERanG, 
WMAP and Planck strongly reassures us of the validity of the so called ACDM cosmological 
model. Nevertheless, we observe that at large angular scales the Universe displays a little 
less anisotropy than we expect [6]. In this paper we analyse quantitatively the corpuscular 
model proposed by Dvali and Gomez in Ref. mi, and show that corrections to the standard 
power spectrnm are larger at large scales, and that the tensor-to-scalar ratio is naturally 
suppressed by cumulative quantum-depletion effects for inflaton masses approaching the 
Planck scale. Althongh conceptnally very different from the standard inflationary scenario, 
the model leads to predictions that satisfy cnrrent experimental bounds. In particular, for 
inflaton masses much smaller than the Planck scale (e.g. around 10^^ GeV), corrections to the 
standard power spectrnm are well within cosmic variance aronnd the standard predictions. 

We will briefly review the original idea of self-gravitating Boson condensates to describe 
black holes in Sec. and show how this model can be applied to inflation in Sec. Ill In 
Sec. W, we nnmerically solve the time-evolution eqnation and show how various cosmological 
parameters behave as functions of time, and how the cosmic microwave background (GMB) 
power-spectrum is affected. Finally, we will comment about our results and outline possible 
fnrther developments in Sec. |V| 

We will mostly use units such that c = 1, h = mpi^, and the Newton constant Gn = 
ip/rrip, where ip and nip are respectively the Planck length and Planck mass, except for 
Section IIVI where we shall also use astronomical nnits. 


^ Not feeling the need to explain the importance of inflationary models, we just mention two of the many 
review papers that have been written on the subject, namely Baumann [T] and Linde [5]. 
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II. THE GRAVITON CONDENSATE 


In recent years, a new way of describing the quantum aspects of black holes was proposed in 
Refs. [9]. The idea is simple and straightforward: black holes can be seen as a self-gravitating 
Bose-Einstein condensates (BECs) of gravitons at a critical point, with Bogoliubov modes 
that become degenerate and nearly gapless, representing the holographic quantum degrees 
of freedom responsible for the black-hole entropy and the information storage [10]. These 
gravitons approximately interact by means of the Newtonian potential 



and correspondingly have an effective mass /i related to their characteristic quantum- 
mechanical size via the Compton/de Broglie wavelength £ ~ h//i = Iptrip/fi. Since they 
are Bosons, they can superpose in an approximately spherical volume of radius i, and total 
energy M = N where N is the total number of constituent gravitons. Within the New¬ 
tonian approximation there is then a value of N for which the whole system becomes a 
black hole. 

This can be seen by considering the effective gravitational coupling constant a = 
~ = /i^/rUp, and estimating the average potential energy per graviton as 

f/~/ xVnW ^-iVa/U • (11.2) 

Black hole components in the depleting region are “marginally bound” when 

+ t/ ~ 0 , (II.3) 

where the kinetic energy is given by Ek — /x. This energy balance yields the “maximal 
packing” condition 

(II.4) 


from which the effective Boson mass fi and total mass M of the black hole are seen to scale 
according to the relations 


/i ~ 


m 


y/N 


= and M = iV /i ~ Vn 


TTlr 


(II.5) 


The horizon’s size i?H and area A are therefore quantized as commonly expected 

Rn ^ \/iV4 




( 11 . 6 ) 


Since this initial intuitive model, other authors have proposed possible developments m 
iniiiMi] and possible cosmological implications [IHl HH] . 


^ One can see that the Newtonian approximation is not far from the general-relativistic behaviour UMM- 
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The Hawking radiation and the negative specific heat spontaneously result from the 
quantum depletion of the condensate for the states satisfying Eq. (II.3). At first order, 
because of reciprocal 2 —)■ 2 scattering inside the condensate, the depletion rate can be 
estimated to be 

1 


^ _ 

s/NL 


(11.7) 


where the first factor comes from the interaction strength a, the second factor is 

combinatoric (there are about N gravitons scattering with other — 1 ~ gravitons) and 
the last factor comes from the characteristic energy of the process. The amount of gravitons 
in the condensate will then decrease according to |H] 


N~-r~--^ + o{N-'). (II.8) 

As explained in Ref. [9], such emission of gravitons reproduces a part of the Hawking radi¬ 
ation which is purely gravitational and contributes to the shrinking of the black hole. The 
usual results are then recovered in the double-scaling limit 77 —)• oo and ip —)■ 0, with h 
kept finite. It is important to remark that, in the present description, particle creation by 
a black hole is not a vacuum process but an actual scattering phenomenon of the back hole 
constituents from the ground state to excited states which display an “effective” Hawking 
temperature. 

Despite its simplicity, the lack of a geometrical description of space-time makes it diffi¬ 
cult to understand the role of any kind of horizon in this model. A possible way to connect 
the BEC approach to the usual geometrical point-of-view of general relativity was later 
proposed [20], using the formalism of Refs. |2T], which introduces the idea of a “fuzzy” 
gravitational radius at the quantum level, starting from the Einstein equation that deter¬ 
mines the classical Misner-Sharp mass for spherically symmetric systems. For example, if 
the quantum state of a single particle is given by a localised wave-packet obtained from 
the superposition of many energy eigenstates, to each eigenstate there will correspond a 
different gravitational radius Rh, with a probability amplitude given by the corresponding 
spectral coefficient. When measuring the particle’s position, we could then observe it with 
some probability inside the expectation value of Rh- Simple evaluations showed that only 
very massive particles, of order of the Planck mass or larger, are likely to appear in a black- 
hole state, but the relative uncertainty in the horizon size is then unacceptably large for 
semiclassical black holes. However, when we likewise analyse a system made of many light 
condensed Bosons IZI, the horizon relative uncertainty decreases with N, which supports 
the idea that large BECs of gravitons at the critical point can be viewed as semiclassical 
black holes in the large-A^ limit^. 


^ If one further considers a scalar field (p coupled to a source J proportional to the field itself, then Eq. (II.61 
is also recovered m- 
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III. QUANTUM INFLATIONARY COSMOLOGY 


As we have already mentioned, the BEC model can be easily extended to cosmology, and in 
this Section we systematically review the scenario introduced in Refs. We shall start 

by noting that, in order to obtain physically sensible results, one hrst needs to reproduce 
the usual background dynamics, as it is described by the Friedmann equation 


2 Stt ip 

H -, 

6 TUp 


(III.l) 


at least to some approximation. We shall in particular show how one can recover separately 
the left-hand side (an effective Hubble parameter) and the right-hand side (the matter 
source) of the above Friedmann equation. The former will be effectively given by the graviton 
condensate, whereas the latter will require introducing an inflaton held 0, which we shall 
further assume satishes the slow-roll condition. Next, we shall show that including purely 
quantum effects leads to a master equation for the time-evolution which will be investigated 


quantitatively in Sec. IV 


A. De Sitter and the Hubble Parameter 


First of all, we note that an observer in the (homogeneous and isotropic) classical de Sitter 
Universe would see a horizon located at the radius 


Re = H-' , 


(III,2) 


where H is the Hubble expansion rate. Like for black holes this radius will determine the 
approximate wavelength and occupation number of the gravitons in the BEC which make 
up the space-time fabric. In fact, let us assume the total gravitational energy conhned inside 
Rh is again E = N fi, where N is the total number of gravitons and fi H/Rn ~ TUpip/Rn 
is the single graviton energy. (Of course, we are here neglecting higher-order corrections to 
this mean held approximation, corresponding to terms that scale with lower powers of N.) 
If we now assume the energy density of this BEC gives rise to an effective curved geometry 
with the Hubble parameter ® 


one must then have 


Sn ip 

^- P 

6 nip 




Sir ip E Stt ip Njj, N 
3 rup R^ 3 nip R^ R^ 


(III.3) 


(III.4) 


Curiously enough, it seems that the radius of the observable Universe coincides with its Schwarzschild 
radius (to a good approximation). 

® For a field-theoretic justification of this step, see Sec. 3 of Ref. [7]. 
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From Eq. (III.2), we again obtain 


Rn ~ y/Nl 


P 5 


(III.6) 


which shows that the critical scalings (II.5) also hold in the cosmological case. The smaller 
the graviton nnmber in the Hubble patch, the stronger the expansion {N oc Also, the 

average graviton density 


p ~ 


M 

K 


rrip 

Nil 


(III.6) 


is qualitatively the same as for a black hole, whose mass is quantised according to Eq. (II.5), 


and where to a larger N corresponds a bigger, less dense and more classical black hole. 

Let us remark that another point in common between black holes and the Friedmann 
models is that both have some affinity with Newtonian gravity. In fact, the full Friedmann 
equation for the scale factor a = a(t) reads 


— 


2GNC'(a) 
3 a 


= k 


( 111 . 7 ) 


where C := dvrpa'^ is constant for dust, and k is the spatial curvature. Eq. (III.7) looks 


exactly like the equation for the conservation of the energy E = 2 pk oi a point particle 
moving along the trajectory r = a{t) in an external Newtonian potential with source mass 
M = C{r). With k = 0 = E for the spatially-flat Universe we seem to live in, this equation 


then yields the same marginally bound condition of Eq. (II.3). It is thus no surprise that 


the de Sitter Universe looks like a black hole turned inside out, where the horizon acts like 
a semipermeable membrane that prevents things from re-entering, whereas in a black hole 
it keeps them from escaping. 

As shown by Dvali and Gomez, there is one important difference between black holes and 
cosmology that we need to take into account now: the passage from de Sitter to Minkowski 
forces us to “introduce a cosmological fluid that plays the role of a homogeneous cosmic clock. 
This role is played by a scalar field, the inflaton <p” [7]. In other words, having recovered the 


left-hand side of Eq. (III.l), we now need to have a suitable matter source in its right-hand 
side. 


B. Slow-Roll Inflation 

In light of the above observation, let us go back to old-fashioned general relativity and 
consider a homogeneous inflaton scalar held of mass m = hXffi and with the simple 
potential 


r(')>) = 2A 


P- 


(III.8) 
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The energy density and the pressure are given by 


P<t> 




and 


= 


- K" 


(III.9) 


respectively. The held (j) evolves according to the Klein-Gordon equation 


3 H d) d) — 0 . 


(III.IO) 


In the slow-roll approximation, |0^| -C |Id| and |0^| 3H(j) ~ the inhaton therefore 

displays a quasi-de Sitter-vacuum equation of state 


P<P - -P<P - “2 • 

This leads to the de Sitter exponential growth of the Universe, with 

An L 




3A^ rup 


The slow-roll parameter 


e : = 


H 

IP 


SiXmHY 3\\ 




= l(^] «1, 


(III.ll) 


(III.12) 


(III.13) 


determines the validity of the approximation, and shows that the Compton length of the 
inhaton must be much larger than the Hubble radius. In particular, we can take the common 
assumption that the exponential inhation turns oh when e approaches 1, or — .Rh- 
In the slow-roll approximation, the total inhaton occupation number inside the Hubble 
patch can be estimated by multiplying the inhaton number density, 

3 A, 


Ptf) ’J 

m 8n 


(HI.14) 

times the Hubble volume ~ 1/H^. By simply substituting in for the previous quantities, 
and recalling Eq. (HI.2), we hnd the total number of inhatons is given by 


AT, ~ H 

~ £2 ^ • 


(HI.15) 


Using Eq. (HI.5), one can then see that the total inhaton number is related to the total 
graviton number via the Hubble constant, that is 


~ ( — A™ i/) iV ~ ^ , 

^ ' Stt y Stt yi ’ 


(HI.16) 


therefore one hnds that 3> iV in the slow-roll regime. 
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Let us now abandon the standard general-relativistic description of cosmology, and return 


to the BEC model of Sec. Ill A We then see that the above relation (III.16) between N and 


can be viewed as the effective background Friedmann equation (III.l), that is 


8 TT f'r 


Sn ir 


wr— 


m 


rur 


(III.17) 


where p ~ rrip/N is the energy density (III.6) of the gravitons in the BEC and ~ 


N^mH^ is precisely the inflation energy density in Eq. (III. 14). Although one must have 
N -C Nfp in the slow-roll regime, the energy densities of gravitons and inflatons remain equal, 
which reminds us that the above effective Friedmann equation is just a statement of total 
energy conservation 

To summarise, we have so far shown how to recover purely classical background dynamics 


compatible with the Friedmann equation (III.l). It is then time to include quantum effects. 


namely the depletion of the background BEC that, for black holes, leads to the Hawking 
evaporation. 

C. Quantum Time Evolution 

Like inside a black hole, we expect that scattering of gravitons will lead to a progressive 
depletion of the background BEC, now further enhanced by the scattering against inflatons 


(see Ref. [7] for more details). If we further assume Eq. (111.17) still holds (to leading order 
in N), given the depletion law for the N gravitons, we obtain the time-evolution of N^, 
which is like implicitly dehning cf), and therefore Ending a way to turn the inflation off. One 
could also introduce a wave-function for the Hubble radius by quantising the Friedmann 


equation (111.17), like one does with the Misner-Sharp mass of a compact source to describe 


the Schwarzschild radius of quantum black holes. The same formalism would yield a “fuzzy” 
cosmological horizon (and therefore H), with an uncertainty that we can anticipate will 
decrease for large N (and, correspondingly, large N^). In fact, we shall apply this formalism 
just to show the effect of excited modes on the background evolution remains negligibly 
small within the slow-roll regime (see Appendix]^. More importantly, we remark that, like 
the depleted quanta make up the Hawking radiation from black holes, the excited inflatons 
will appear to an observer as the scalar density perturbations, and the excited gravitons as 
the tensor perturbations. Like for black holes, in this model there is no production of quanta 
from the vacuum, and the fact that the higher energy of the quanta above the ground state 
does not significantly affect the background dynamics can be translated into the statement 
that the backreaction of density and tensor perturbations is mostly negligible. 

In order to study the time evolution of our inflationary background, we can directly 


employ the relation (HL3) between the Hubble rate H and the energy density p of the 


® In general relativity, Eq. (III.l) is the Hamiltonian constraint that follows from the Einstein-Hilbert action. 













gravitons. The novelty with respect to Section IIIB , where the inflaton was introduced in 
order to recover the slow-roll evolution, is that we also want to take into proper account the 
effect due to the depletion of the BEC. In particular, the master equation for the number 
N of gravitons in the BEC used in Ref. [7j can be obtained rather straightforwardly. 

The rate of change of N will in fact contain two contributions: a “classical” (mean field) 
one due to the presence of the inflaton and described by the (effective) Friedmann equation, 
and a purely “quantum” one due to the depletion of the background BEC, 


N 

N 


Ncl Ndep 

^ ~w~ 


(III.18) 


The former contribution is obtained via the Friedmann equation (III.17) and Eq. (III.6), 
that is 


Nd 

N 


P 

P 


1 dTT^ ^ H _ 

Ip ~6r ~‘^H~ 


-2eH . 


(III.19) 


The quantum depletion is again given by iV^ep ~ —T, where T is obtained from the expression 
in Eq. (II. 8[) with a small change to account for the presence of more species. As explained 


in Ref. we must consider three kinds of interactions: graviton-graviton {g — g), graviton- 
inffaton {g — (p) and inffaton-inffaton (0 — 0). Since the inflaton wavelength is virtually 
infinite in the slow-roll approximation (A^ 3> Rr), the 0 — 0 process is irrelevant with 
respect to the scatterings involving gravitons. The g — g and g — (p scatterings differ by a 


combinatoric factor which is by far larger than N in the slow-roll regime [see Eq. (III. 16) 

It is therefore the g — 4> process that dominates and depletes gravitons and inffatons at 
the same rate T^^. Taking furthermore into account that the characteristic momentum 
transfer is h ih ~ £p rup H, we find that the depletion rate Tg^ can be well approximated by 
(c.f. Eq. (5.12) in Ref. [7]) 


Tg^c^a^N^NH . 


(III.20) 


Here, the first factor (the square of the dimensionless self-coupling a = {IpHY = (f'p/Rn)^) 
stems from the fact that the dominant process is again give by 2 —)■ 2 scatterings. The 
factor Nfj, N counts the available pairs in the g — p scattering, and H is the characteristic 
momentum transfer, as mentioned above. Using the “maximal packing” condition ~ 1, 
we have 


7V = iV^~-r,^~-i7«2iV^iV 


With Eq. (III. 16), this leads to 


77dep Tg^ 1 1 

TWIpl^N 


1 iV0 
y/Nlp N 


73 H 


(III.21) 


(III.22) 
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Putting together these pieces, and temporarily omitting numerical factors in order to high¬ 
light the functional dependence on N, we hnally hnd the master equation (c.f. Eq. (5.27) in 
Ref. [7j) 


N 

N 


~ e if — 


After replacing e ~ m'^/{3 and H 

dependence manifest, we get 


H 

7 ^ ■ 

- rtjj _ 


(III.23) 

{'/N ip)“^, and making the time- 


• ^ 2 _ mp 

3 ’^p f'p ipmN (t) 


(III.24) 


which reproduces the classical Friedmann dynamics in the double-scaling limit N ^ oo and 
h —)■ 0, with Gn and m kept hnite. 

Corrections arise when going to higher orders in 1/iV, which hnd correspondence in higher 
orders in the BBGKY-hierachy |22]- The inhuence of the hrst of such terms, i.e. the depletion 


term in Eqs. (III.23) and (III.24), gives rise to quantum corrections, which differ from the 
standard semi-classical results to which we will compare our respective hndings below. All 
these corrections cumulate, and higher-order terms will have to be taken into account after 
a certain time. This is precisely what limits the total duration of inhation (c.f. Sec. 5.4 
of Ref. m) as the mean-held (background) description becomes poorer and poorer. Then, 
eventually, one would need to go beyond the hrst 1/N term, or, respectively, the second-order 
truncation of the BBGKY-hierachy. 

We again remark that we have simply retained the functional dependence on m and 
N, and set a numerical coefficient in front of the second term on right-hand side equal of 
Eq. (III.24) to one. This simplihcation should be kept in mind when discussing the results in 
the next Section, since a diherent relative numerical coefficient between the two contributions 
to the evolution of N will likely change the hnal hgures, but this analysis is left for future 
investigations. 


IV. COSMOLOGICAL PARAMETERS 


Of course, Eq. (III.24) has to be solved with appropriate boundary conditions, which we will 


choose to specify at the end of inflation t = tend, i-e. N(tend) = W*, for a certain number of 
condensed background gravitons A*. This will be estimated, given the energy scale m, by 
{mp/m)‘^] for instance for m = 10^^ GeV one has A* ~ 10^^. We then solve Eq. (III.24) 




numerically and backwards from tend over a time interval [tstart, tend] corresponding to a total 


number Ae = 200 of e-foldings. For the simple inflaton potential (III.8) one has 

tend)^ , 

where the time tend is characterised by the condition e(tend) = 1- 


m 

Ae = (tstart 

0 


(IV. 1) 
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The upper panel of Fig. shows the time-dependence of the number of condensed gravi¬ 
tons for various values of the inflaton mass, and we display the relative correction of the full 
solution to Eq. ( III.24[ ) with respect to the one without quantum depletion in the bottom 
panel of the same Fig. We observe that the shapes for the different values of the mass 
are very similar, with the quantum-depletion correction becoming increasingly important 
for increasing mass. In fact, this correction roughly scales as m^. ^ 

The top panel of Fig. [^displays the slow-roll parameter e(W) as a function of the number 
of e-foldings N^. It is easy to see that e is larger for larger mass m at earlier times, and it 
then approaches one at the end of inflation (which corresponds to W = 200 here). Simi¬ 
larly to Fig. [T| we plot in the bottom panel of Fig. the corresponding relative correction 
^corrected/^uncorrected ” 1, which is also approximately proportional to m^. The same holds true 
for the relative correction hfcorrectedZ-ffuncorrected — 1 of the Hubble parameter (which scales 
like 1/v^), being presented in Fig. 

At this point, one might wonder about the corrections to the tensor-to-scalar ratio, since it 
is strongly sensitive to the inflaton mass, but is also expected to undergo quantum-depletion 
corrections. This is even more true as it probes the total duration of the inflationary expan¬ 
sion, and is approximately given by (c.f. Eq. (5.36) of Ref. [7]) 




1 H{t) 

Atart 


(IV.2) 


For evaluating this quantity we employ results from the Planck paper [23] ; we take the pivot 
scale fcpivot = 0.05Mpc“^, and compute r(tpivot) at the time tpivot when /cpivot exited the 
horizon. Of course, as the Hubble parameter is affected by quantum depletion, we remark 
that tpivot also acquires a correction. 

For the potential used in this work, the Planck collaboration gives the upper bound 
’" 0.05 ^ 0.12 (Planck-I-WP 95% CL, c.f. Tab. 9 in Ref. [23]) More recently, this bound 
was confirmed after including BICEP2 and Keck Array data in Ref. [21]. We tabulate the 
specihc values of r at the chosen pivot scale in Tab. One can observe that masses up to 
m ~ 10“^mp are consistent with this Planck bound. Furthermore, we see that, for large 
masses, the influence of the quantum depletion is rather strong and it actually lowers the 
value of r, hence, it generically broadens the range of compatible masses. This effect is 


^ Of course, for large-enough masses and for a sufficiently long duration of inflation, the underlying approxi¬ 
mations will always break down at some point, as quantum-depletion effects continuously and cumulatively 
lower the amount of condensed quanta, making successively the l/7V-corrections more and more impor¬ 
tant. The larger the mass the sooner this will happen. However, for all the cases studied in this work, 
those corrections are under control. 

® The bound obtained by including the running of the spectral index is r < 0.25 (ACDM + r + drig/d ln(A:), 
Planck-hWP, c.f. Tab. 5 in Ref. [23]) 
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Figure 1: Upper panel: The number of condensed gravitons N{Ne) as a function of the number of 
e-foldings for various values of the inflaton mass m = 10 “^’rup (bottom to top). Lower panel: 
The relative correction iVcorrected(.^e)/.^uncorrected(-^e) — 1 of the quantum-corrected to classical 
evolution of the number of condensed gravitons as a function of the number of e-foldings, for 
various values of the inflaton mass m = 10“^’■■■’rup (top to bottom). Masses are in units of nip. 
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m[mp] 

ratio 

10"® 


10"^ 

10"3 

5•10-3 

10-2 

1.5 • 10-2 

2 X 10-2 

corrected / uncorrected 

1 

1 

1 

0.93 

0.43 

0.25 

0.18 

0.14 


Table I: The tensor-to-scalar ratio r at the pivot scale fepivot = 0.05 Mpc red numbers indicate 
violations of the bound r < 0.12 (see main text). 

stronger the higher the inflaton mass, as expected since the effect of quantum depletion 
increases with increasing mass (c.f. the above discussion). 

Let us next come back to the quantum-depletion effect on the Hubble parameter. As can 
be observed in Fig. H will be suppressed more strongly at earlier times. This immediately 
leads us to investigate whether we can learn something about the observed suppression of 
the low multipoles [6] of the CMB temperature auto-correlation function. 

In general, the spectrum of (relative) temperature fluctuations 6T/T can be described by 
an infinite set of correlation functions 

C(9i2. 9iS. . . . := %(k) ■ ■ ■ %(ln)) (IV.3) 

with n G N and the angle brackets denote an average over all possible directions li, for a given 
configuration of angles 9ik '■= arccos(Zi ■ Ik) (see Ref. [22] for an extended review). Assuming 
approximate Gaussianity, the spectrum can be described by the two-point function 

C(«):=C(e,2) = ^T((j)^(y^^ (IV.4) 

alone. It is convenient to decompose it as 

2 / -I- 1 

C(9) = C, P,[cos(9)], (IV.5) 

1=2 

where the multipole moments G IR are expansion coefficients, is a Legendre polynomial 
of degree £, and we subtracted the monopole and the dipole. Again, under the assumption 
of homogeneity and isotropy, the average over directions — appearing in the definition of 
Off) — is equivalent to the average over different points, keeping the directions li and I 2 
fixed with the condition 6 = arccos(Zi ■ I 2 ), for a given angle 9. 

The multipole moments are, modulo transfer functions and prefactors, given by a 
fc-integral over the inflationary power spectrum which is approximately proportional to 
(c.f., e.g.. Ref. [I]). To this end, in Fig. we plot the quantity 

M'S /L') 

A(fc) := 1 - J ^ 

^ uncorrected v ^ / 
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Figure 2: Upper panel: The slow-roll parameter e(A^e) as a function of the number of e-foldings 
for various values of the inflaton mass m = 10“^’(top to bottom). Lower panel: The 
relative correction ecorrected/cuncorrected — 1 of the quantum corrected to the classical evolution of 
the slow-roll parameter e{Ne) as a function of the number of e-foldings, for various values of the 
inflaton mass m = 10“^’■■■’rup (top to bottom). Masses are in units of nip. 
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Figure 3: Upper panel: The Hubble parameter H{Ne) as a function of the number of e-foldings for 
various values of the inflaton mass m = (top to bottom). Lower panel: The relative 

correction Hcorrected/Huncorrected — 1 as a function of the number of e-foldings, for various values of 
the inflaton mass m = 10“^’rup (top to bottom). Masses are in units of nip. 
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Figure 4: The quantity A{k) := 1 — -f^corrected(^)/-^uncorrected(^) ^ function of the wave-number 
k for various values of the inflaton mass m = 10“^’“^’“® rup (top to bottom). Masses are in units 
of nip. 


which shows us exactly how the correction to the integrand of the mentioned C^-integral 
behaves. For small scales it is approximately proportional to m^. 

For the calculation of the coefficients we use a modified version of GAME |2S], employing 
the CMBQuick package ^7\. Fig. shows the ratio of the uncorrected to the quantum- 
depletion-modified multipole moments. The important result here is given by the green 
dashed line corresponding to m = lO^rUp, which, for the potential (III.8), is the standard 
required value leading to the normalisation of the matter power-spectrum used by the Planck 
collaboration |23] (c.f. p. 64 in Ref. pQ). For illustrative purposes only, we also display the 
relative corrections arising from the choice of m = 2.5 ■ 10“^ mp, which we add to the 
standard (properly normalised) result. Although the base curve onto which the corrections 
are added would be different, the resulting line illustrates the behaviour of the respective 
quantum-depletion correction. We indeed observe a potentially significant suppression of 
the multipole moments for small £, which corresponds to large angular scales. As expected, 
this suppression is stronger for larger values of the inflaton mass. 

In Fig. we present the full power spectrum £{£+!) C(,/{2ti) Tq , where Tq is the CMB 
temperature today. Here, this suppression is visible for small £, albeit it would only be 
notable for rather large values of the inflaton mass. To this end, and for illustrative purposes 
only, we again show the relative change to the power spectrum for m = 2.5-10“^ rup, bearing 
in mind that the normalisation for this case would not match the Planck one. Note that. 
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Figure 5: The ratio of the uncorrected to the quantum-depletion modified 

multipole moments. Masses are in units of nip. 

even for this rather high value, all corrections remain well within the cosmic-variance band. 
Of course, the precise shape of this band is slightly model dependent. 

Let us investigate whether this effect might also tell us something regarding the seeming 
mismatch of the measured to the standard ACDM-calculated angular temperature auto¬ 
correlation function C(9) (c.f. Ref. [28] and references therein). Therefore, in Fig. we 
compare C{6), which is computed by summing multipole moments up to £ = 800, for var¬ 
ious cases (with and without quantum-depletion modihcations) to the Planck data. We 
observe that the larger the mass, the more the corrected curves approach the one inferred 
by the Planck measurements. Although this effect seems too small to accommodate for the 
aforementioned difference between measurements and the standard theory, we think that it 
deserves further investigations. 
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V. DISCUSSION AND OUTLOOK 

In this work, we have quantitatively investigated the full time evolution of inflationary space- 
time in the corpuscular description proposed by Dvali and Gomez [7]. We have computed 
several cosmological quantities, like the slow-roll parameter e, the Hubble parameter H, and 
the tensor-to-scalar ratio r. For all those quantities, we quantihed the intrinsic quantum 
effects inherent to the studied corpuscular structure of gravity. These effects are actually 
stronger the longer inflation lasts, and increase with increasing inflaton mass. 
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Figure 6: Comparison of the power spectra £(£+l)C£/(27r)TQ for the standard case (black, solid), as 
well as for the quantum-depletion corrected results (m = 10“®mp, green, dashed; m = 2.5-10“^mp, 
red, dot-dashed). The shaded region depicts the cosmic-variance band. Masses are in units of mp, 
and the case m = 2.5 • 10“^mp is included for illustrative purposes only. 


Regarding the tensor-to-scalar ratio r, we found that the quantum-depletion effects could 
be relevant, as r is sensitive to the total duration of the inflationary phase and these effects 
are cumulative in nature. In particular, r would be lower than in the standard semi-classical 
computation for inflaton masses approaching the Planck scale. 

We found that quantum depletion decreases the Hubble parameter, whose amount we 
precisely quantihed for different inflaton masses. The gravitational potential is, in turn, also 
suppressed, and so are the CMB multipole moments. We computed the precise amount of 
this suppression for several values of the inflaton mass, and found that it is strongest for 
low multiples, or equivalently, on large scales. The suppression is also stronger again for 
increasing inflaton mass, and is completely consistent with observations, in particular for 
the CMB multipoles and the temperature auto-correlation function. 

Having studied the simplest inflationary model, consisting of a single real scalar held with 
a canonical kinetic term, as well as a mass term, we are tempted to numerically investigate 
corpuscular effects on more rehned models of inhation, although we hnd it remarkable that 
already the simple set-up studied here, where just the functional dependences in the master 
equation (HI.24) are retained, is capable of reproducing the standard results. 

On the qualitative level, we believe that our hndings are generic, as the occurrence of the 
“quantum clock”, originating from the quantum depletion (and working against the semi- 
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Figure 7: Comparison of the angular temperature auto-correlation function C{9) for the standard 
case (black, solid), the quantum-depletion corrected results (m = 10“®mp, green, dashed; m = 
2.5 • 10“^mp, red, dot-dashed), as well as for the Planck data (blue, thick, dotted). Masses are in 
units of mp, and the case m = 2.5 • 10“^ mp is included for illustrative purposes only. 

classical one) is a necessity in the corpuscular picture (cf. Ref. |8]), and the calculation of 
the relevant 5-niatrix elements do not really depend on the specihc choice of the inflaton 
potential (cf. Ref. [21] )• The direction in which the corrections are going should always be 
the same, and they should be stronger for lower multipoles (which correspond to earlier 
times of horizon crossing). However, on the quantitative level, different potentials, like the 
monomial ones considered by the Planck collaboration 1211 , will lead to different strengths of 
the corrections. On a related note, but investigating the question of the occurrence of eternal 
inflation in the corpuscular picture, it was shown in Ref. [30] that monomial potentials lead 
to qualitatively similar effects, and these will be stronger the steeper the potential. 

Finally, it would be very interesting to further investigate the cumulativeness of the 
quantum corpuscular effects, which seems likely to both tighten and loosen constrains on 
inflationary and possibly modihed- / quantum-gravity models. 
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Appendix A: Effect of Non-Condensed Quanta 


The effect of condensed graviton quanta is to generate an effective Hubble expansion pa¬ 
rameter H ~ that satishes the usual Friedmann equation. Their occupation number 

N is partly driven by the presence of the inflaton, and by the depletion, which result in 


the master equation (III.24). However, if we want to take into account the effect of the 


non-condensed quanta, we should also add their energy density to the energy density of the 
“ground state”, like it was done for black holes in Ref. EOI- In fact, the depletion causes 
the ground-state energy to decrease by reducing N, but the quanta scattered off the ground 
state have higher energy, which should affect the Hubble radius. 

Let us then replace p with p -|- 5p, where 5p is the energy density of depleted gravitons. 
This will of course affect the Hubble parameter ~ p and the slow-roll parameter e ~ 
rn ?The perturbation 5p is already provided in Ref. [7], and is given by 


5p = 


6Nx mj. 


p'-p 


R3 R 


(A.l) 


H 


where Nx is the number of deplete particles of wavelength A, both gravitons and inflatons. 


From Eq. (HI.21) we can see that 6N^ = 6N = N^/N, which gives us 


^p = ra i — ^ 


(A.2) 


Since the slow-roll parameter is also expressed via e = —iL/iP^, we hnally obtain the dy¬ 
namical equation 




rrir 


TUr 


+ 6mN^/‘^{t)/mp) ipUiNit) 


(A.3) 


Note that Eq. (A.3) precisely reduces to Eq. (HI.24) for large N. 

As before, we shall use the quantity A(P) := /Lfuncorrected(^) ^ measure of 

the relative strength of a given correction. In Fig. we then display the ratio Afuu/Anwroniy, 


where Afuu is computed from the evolution given by the full Eq. (A.3), whereas Anwroniy 
is obtained by subtracting the quantum-depletion term from that equation. ® We observe 
that this ratio is essentially independent of the chosen inflaton mass and is always much 


® AHWFoniy thus represents the correction due to the higher energy of the excited quanta only, which would 
be obtained by applying the formalism of the horizon wave function [21] to the present case. 
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Figure 8: The ratio Afuu/AHWFoniy as a function of wave-number fe; masses in the Figure are in 
units of rup. 


larger than one. This shows that the correction omitted in the analysis of Sec. |lV]is indeed 
negligible, and one can jnst employ the master eqnation 


2 m? rrip 

3 TTip £p m N{t) 


(A.4) 


In Fig. 1^ we illustrate the influence of the various quantum corrections on the time- 
dependence of if, for various quantum corrections. 

Let us hnally note that the de Sitter condition FT = 0, or, = 0 can be obtained for a 
number of graviton quanta in the BEC of order 


N ~ 



(A.5) 


which, for the typical inflaton mass m ~ 10^^ GeV, would be of order 10^. 
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